> 

m 

cn 

o 

as 



X 



BUHEP-99-25 
McGill-99/31 
OITS-681 
SHEP-99-16 



On the QCD Ground State at High Density 

"cK 
a^ 

^ I Nick Evans^, James Horinuzdiar^, Stephen D.H. Hsu^, Myck Schwetz"^ 

Q , ^Department of Physics, University of Southampton, Southampton, S017 IBJ, UK. 

N—/ . "^Department of Physics, McGill University, 3600 University, Montreal, Quebec, H3A 2T8, USA. 

^_j . '^Department of Physics, University of Oregon, Eugene OR 97403, USA 

Department of Physics, Boston University, Boston, MA 02215, USA. 



October, 1999 



Abstract 



Ph! 

D , We investigate the possible ground states of QCD at asymptotic densities, where the 

'■^ theory is expected to exhibit color superconductivity. We characterize the color-flavor 

structure of possible diquark condensates, and find those that are energy extrema by 
solving the weak-coupling Dyson-Schwinger equations, including Landau damping and 

5r , the Meissner effect. We show that, as previously anticipated, in the two fiavor case 

the vacuum breaks SU(3) color to SU(2) and in the three flavor case the vacua with 
color-flavor locking (CFL) have the lowest energy. We identify a number of relatively 
fiat directions in the potential along which the pattern of gauge symmetry breaking 
changes and parity is violated. We discuss possible phenomenological consequences of 
our results. 



1 Introduction 

QCD at high density and low temperature is a color superconductor |[I|-[^ characterized by 
the formation of a diquark condensate in the attractive 3 color channel. The condensation 
is analogous to Cooper pairing in ordinary superfluids or superconductors, and occurs even 
via arbitrarily weak attractive interactions due to the presence of a Fermi surface [^, ^, ^ 



Recently it was discovered [^ ^, |T^, |T8], ^ that long range magnetic fluctuations enhance 



the condensation, leading to a gap which behaves as 

A~/xf7-5exp(--^) (1) 

in the weak coupling (small g, or large fi) limit. In this limit (|T|) is the gauge invariant, leading 
order result of a systematic expansion in powers of g. The properties of the condensate are 
easy to determine in the case of two quark flavors. Because the condensate occurs between 
pairs of either left (LL) or right (RR) handed quarks in the J=L=S=0 channel |^, and the 
3 color channel is antisymmetric, the quarks must pair in the isospin singlet (ud - du) flavor 
channel. However, even in this case there is a subtlety, as the relative color orientations 
of the LL and RR condensates are not determined by the usual leading order analysis. A 
misalignment of these condensates violates parity, and further breaks the gauge group beyond 
SU(3)c —>■ SU(2)c. As we will discuss below, an analysis of the Meissner effect is necessary to 
determine the relative orientation. There are thus a number of unstable configurations of only 
slightly higher energy with different color-flavor orientations (and hence different symmetry 
breaking patterns), leading to the possibility of disorienting the diquark condensate. We 
include a discussion of possible phenomenological signals associated with these phenomena. 

The generalization to three flavors is far from straightforward. Again, one can show 
that the condensate must occur in the J=L=S=0 and color 3 channel. The Pauli principle 
then requires that the flavor structure again be antisymmetric ~ {qiqj — qjQi), for quarks of 
flavor i,j. Thus, one can have combinations of condensates which are in the 3 of both color 
and flavor SU(3)l or SU(3)r. Due to the chirality preserving nature of perturbative gluon 
exchange, there is no mixing of LL and RR condensates, which form independently. One can 
immediately see that there are a number of possibilities. For example, the condensates for 
the three flavors and both chiralities might all align in color space, leading to an SU(3)c -^ 
SU(2)c breaking pattern. A more complicated condensate has been proposed [|1^, called 
Color Flavor Locking (CFL), in which the 3 color orientations are "locked" to the 3 flavor 
orientation. 

In this paper we determine the nature of the energy surface governing the various color- 
flavor orientations of the condensate. Let us begin by characterizing the color-flavor conflg- 



uration space of condensates. We consider the ansatz 

Aab L,R AC L,R ^abc^ /oN 

ij = ^k e eijk , [Z) 

where a,b are color and i,j flavor indices. L and R denote pairing between pairs of left and 
right handed quarks, respectively. Under color and flavor A transforms as 

A^ -^ U.A^V^ , (3) 

where Uc is an element of SU(3)c and V^ of SU(3)l. A similar equation holds for A^. It is 
always possible to diagonalize A^ by appropriate choice of Uc and V^: 

f a 0\ 



A' 



(4) 
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V c 7 

Generically, there does not exist a V^ which diagonalizes A^ in this basis. In the CFL 
case, where the diagonalized A^ is proportional to the identity, a = 6 = c, it is easy to 
show that one can choose V^ such that A^ = ±A^. These two configurations are related 
by a U{1)a rotation (see section 3). Hence, they are degenerate in the high density limit 
where gluon exchange dominates. Instanton effects, important at intermediate density, favor 
A^ = A^. Note that parity, if unbroken, requires A^ = A^, and hence implies simultaneous 
diagonalizability. 

In what follows we consider the potential vacua parametrized by a,b,c. First, we use 
the Dyson-Schwinger (gap) equation to determine which of these configurations are energy 
extrema. Next, we compute the energies of the extrema to determine the true groundstate. 



A similar analysis has been carried out by Schafer and Wilczek [|I6| in the approximation 
where gluon interactions are replaced by local four fermion interactions. They concluded 
that the CFL vacuum had the lowest energy. Here, we include the gluons in the analysis, 
introducing long range color-magnetic fluctuations (controlled by Landau damping) and 
Meissner screening into the gap equation and vacuum energy calculations. 

We find that the CFL vacua remains the lowest energy state, at least at asymptotically 
high densities where the calculation is reliable. The Meissner effect is a small correction to 
the vacuum energy at asymptotic densities. At lower densities where the gauge coupling 
is large the Meissner terms become more important and tend to disfavor CFL relative to 
the absence of a condensate. We do not know whether they are ever sufficient to remove 
the superconducting phase but they will lower the energy difference between the vacuum 
and unstable condensates with different color and flavor breaking patterns. Configurations 
which satisfy the gap equations but are not the global minimum of energy are presumably 
saddlepoints, since they are continuously connected to the CFL vacuum via color and fiavor 
rotations. 



2 Gap Equation 

In this section we determine the subset of parameter space for which our ansatz (^) satisfies 
the gap equation. Because the gap equation results from the extremization of the effective 
action, its solutions are energy extrema. At asymptoticaly high densities (weak coupling) the 
diagrams (a)-(c) in figure 1 give the leading approximation to the effective action. Note that 
in these diagrams the quark propagators include the diquark condensate (see ( p!OD below), 
and the gluon propagators include Landau damping, but not the Meissner effect. The latter 
arises from the condensate-dependence of quark loops in diagrams (c) and (d). The resulting 
gap equation (figure 2, with condensate shown exphcitly at lowest order in A) is given by 

S~\q) - So\q) = ^9'J-0^, r^ S{k) rf DZik - q) , (5) 

where 




r^ = ^^ (6) 





D^Q is the gluon propagator, include the effects of Landau damping and Debye screening 
(we assume Feynman gauge throughout): 

where Pt,l are transverse and longitudinal projectors. The analytic forms of F and G are 



given below in (^). The small — expansion of G leads to the Landau damped magnetic 



gluon propagator 

pfiu 

Dli^iqo.q) = ^ r^ , (8) 

while the expansion of F leads to the usual longitudinal propagator, with Debye screening: 

'"D — ^^/ 27r2 • 

We will restrict the color group structure in the gap equation to the attractive anti- 
symmetric 3 channel: 

^aM - liSaAd-SatSed) , (9) 

which projects out the anti-symmetric part of S{k) in color space in the gap equation. Here 
S is the fermion propagator for the spinor (7/'^, ip^^) with i a flavor index and a a color index. 

For the three flavor case S can be written explicitly as an 18 x 18 matrix in color flavor 
space. The inverse propagator may be written 







Figure 1: Vacuum energy diagrams 





Figure 2: Dyson- Schwinger equation 
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Because we are deahng with a diquark condensate the non-trivial part of the gap equation 
involves the lower left 9x9 block. We will refer to this sub-block of the propagator S as 5*21. 

For a particular ansatz A to be a solution to the gap equation we require that the 
color antisymmetric part of T"^ S2i{k) T^ (corresponding to the 3 channel) be proportional 
in color- flavor space to S^^{q) — S^'^^q) = A(g), which appears on the LHS of the gap 
equation. This requires some justification, as the matrices that appear on the RHS of the 
gap equation appear inside the integral. If, for simplicity, we restrict ourselves to ansatze 
which correspond to constant color-flavor matrices times a function of momenta, then this 
condition is implied. In principle, there can be more exotic solutions in which color and 
flavor orientations rotate in momentum space, however it seems unlikely that such solutions 
exist. We note that the equality must hold for all values of the external momentum q, and 
that the set of functions D{k — q) are likely to form a complete basis for functions of k, since 



they are essentially smeared delta functions of (k-q). Thus by taking appropriate linear 
combinations of the gap equation we can see that T"^ 5*21 (A;) T^ must be proportional to 
A when integrated against any arbitary function of k. Hence the proportionality must hold 
without the integral. 

The propagator may be found by inverting the sparse matrix in ( ]T0| ) using Mathematica. 
Only three ansatze satisfy our condition: a = b = c; a = b,c = 0; b = c = 0. We refer 
to these solutions as (111) (color-flavor locking), (110) (3^0 breaking) and (100) (3^2 
breaking) respectively. 

For these ansatze the color antisymmetric part of T'^S2i{k)T^ has the form of a constant 
multiplying the matrix form (|TT|) with a, b, c set to or 1 as is appropriate for the ansatz. 
The constants are (here /^ = (|A;| — /x)^): 

2A {kl -f + 3A2) 



;iii) 



(fc2-/2 + A2) (fc2-/2 + 4A 



2^ 



A A 



(A;2_/2 + A2) (A;2-/2 + 2A2) 

The integral over / can be performed by contour integration, yielding the following gap 
kernels 



, , 2 A 1 A 

(111) : - , „ . +- 



;iio) 



3 ^kl + A2 3 ,y^2 + 4A2 

A A 



2jkl + A2 2Jkl + 2A2 



(100) : ^=^= (13) 



'kl + A2 
Let us now simplify the gap equations. We neglect go, as compared to |g|, as well as 



anti-particle contributions (see, e.g., |]T8[, for details). We obtain 



A(Po) = -^ fdqofd cose (- , , i ,^'"'^ ,,,,,, ,, (14) 



127r2y ^''J yi-cose + {G + {po-qoy)/{2fi^) 



i + I cos \ J.., . 



where 



l-cose+{F+{po-qo)^)/{2^i^), 



n 2?Vt iQo^ (m\\ ^,. 1, (x+l 



G 



■m 



2 m 

k1 




(15) 



and K{qQ) is one of the gap kernels from 



These gap equations can be solved numerically. We first present solutions neglecting 
the Meissner effect. The results for A vs po are displayed in figures 3 and 4 for the three 
ansatze[| (The spatial momentum p is taken to lie on the Fermi surface.) The curves lie very 
close to each other but as we will see below give quite different contributions to the vacuum 
energy. Note that the gap solutions we obtain have broad support, from the Fermi surface 
to l,kQ ~ yU. However, this is likely a consequence of the approximations used in (0), in 
which all momenta are assumed to lie close to the Fermi surface. 
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Figure 3: Gap Solutions for \i = 400 MeV 



A complete analysis must also include the Meissner effect, that is the screening of the 
gluons induced by the formation of the gap. The leading order contribution to the gluon 

, although the shapes of the curves are 



^Our results differ somewhat in normalization from those of 
in agreement. 
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Figure 4: Gap Solutions for /x = 10 MeV 

vacuum polarization P{ko, k) conies from the off-diagonal terms in the fermion propagator- 
matrix 



5Ptiu{qo,q) = g" 



(M 



rTr 



7,Su{k)CY.C-'S2i{k-q) 



(16) 



where we have dropped the group theory factors. It is convenient to use the formalism of 
1^ which exploits simplifications due to decomposition of the fermion propagator into a sum 
of projections onto different chirality and helicity channels. After some tedious algebra one 
gets for (O) 



SPi,iqo,q) = 2g' 



d^k |Ap gij + ^{kik - Qj + kjk - q,- 



(17) 



(27r)4 ki + e{ky {ko - qoY + e{k - qf ' 

where k = k/\k\ and e(A;)^ = (|A;| — /i)^ + |Ap. It is easier to compute the contribution to 
the magnetic gluon mass G directly applying the transverse projector P^ = (5ij — qiqA to 
the gluon vacuum polarization while using the HDL approximation (the momentum in the 



loop k r^ fi and is much bigger than the momentum transfer q) |2 

SG{qo,q) = -P[jPij{qo,q) = / / TTTZU ,. /,.N2 71 „ no , .,. „A2 • (l^) 



Further simphfication comes after switching on a small temperature and performing a sum- 
mation over frequencies |^. Note that, because the system is already decomposed into 



particles and anti-particles about the Fermi surface, one should apply the summation for- 
mulae as if // = 0. Finally, one finds 

<^^[qo, q) - ^2vr)3 J e{k) e{k - q) ql + (e(fc) + e{k - q)Y " ^''^ 

One can see, either by analytical approximation or numerical evaluation, that 6G{qo, q) is of 
order r7i|, for go ~ Q' ~ ^) and falls off like 1/go or 1/g as either become large [|r^. While 
this is of the same order as Landau damping, numerical evaluation shows that the Meissner 
contribution is somewhat smaller. 

As we are only interested in the size of the contribution of the Meissner effect, we use 
the following approximation, which is an overestimate of the effect: 

(5G(go,g)^m^ , ^l ^ , (20) 



q^ + ql + A, 



where Aq is the maximum value of the function A(/i;o, k). The gap equations were numerically 
solved for all three gap kernels, and the results are shown in figures ^ and |^. The effect is to 
decrease the size of the condensate but it is a small perturbation on the solutions obtained 
previously. 



3 Vacuum Energies 

To determine which of the above gaps is the true minimum energy state we must calculate the 
vacuum energy, which receives contributions from vacuum to vacuum loops of both quarks 
and gluons (figure 1). We use the CJT effective potential, which is a function of condensates 



Hi: 



74 

V{S,D) = -^|^[trln5(p)/5o(p) + tr(l - 5(p)/5o(p)) - ^tr lnD(p)/Do(p) 

- '-tT{l-D{p)/Do{p))] + IJ J ^,-0^,[^'^S{p)TS{p+k)D{k) 

+ ■■■ (21) 
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Figure 5: Gap Solutions for fj, = 400 MeV 

where for convenience we suppress appropriate color, flavor and Dirac indices. 5*0 and 5* 
correspond to bare and full fermion propagators, Do and D to bare and full gluon propagators 
and r to full vertices. The ellipsis denote gluon self-interaction loops and terms which are 
higher order in g. In our approximation, which is essentially Hartree-Fock (lowest order in 
coupling), the F's become bare vertices. 

Extremizing with respect to appropriate propagators and vertices one obtains a set of 
gap equations. The fermion gap equation is the one we studied in the previous section, while 
the gluon gap equation produces Landau damping. We wish to compare values of V{S, D) 
corresponding to our three solutions to determine which one is the true vacuum^. It is easy 
to show that the value of the effective potential evaluated on the gap solution {S*,D*) in 
the Hartree-Fock approximation is given by: 



V{S\D*) 



(2^ 



tihiS{p)/So{p) 



(22) 



^The difference in energies V will be gauge invariant, whereas actual values are not. 
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Figure 6: Gap Solutions for /x = 10 MeV 

Diagramatically, this is equivalent to the graph of figure 1(a) when evaluated on the gap 
solution. 

The fermion loops are most easily calculated by going to a basis where SqS^^ is diagonal 
in color-flavor space. Note that the gap matrix A has non-trivial Dirac structure that must 
be accounted for |^: 

A = A175P+ + A2^,P. , (23) 

where P± are particle and anti-particle projectors. Our analysis has been restricted to the 
particle gap function Ai. The anti-particle gap function A2 has its support near ko ~ 2/i, 
and its contribution to the vacuum energy is suppressed. There are 18 eigenvalues, which 
occur in 9 pairs. The product of each pair is of the form 

A\ko,k) \ 



- [l + a 

\ A;g + (|fc|-/i)V 

where a is an integer. For our three cases we obtain the following sets of eigenvalues: 

(111) -^ 8 X {a = 1} , 1 X {a = 4} 



(24) 



11 



(110) - 


4 X {a 


= 1} , 2 X {a = 2} 


(100) ^ 


4 X {a = 1} 


The binding energy is of order 




E, ~ -/ 


i^kdko In 


[l + a ^'^'-^) 1 


~ —a /J 


2a2 
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(25) 



(26) 

(27) 

where Aq is the maximum value of the gap function A(A;o, /c), which has rather broad support 
in both energy and momentum space away from the Fermi surface, extending to kQ,k ^ fi. 
A more precise answer than (|2^) requires numerical evaluation, but it is clear that the result 
scales with a and has only a weak (logarithmic) dependence on the variations in the shape 
of A(/co, k). Substituting our numerical results for the gaps in the three cases, it is easy to 
establish that 

E(lll) < E(llO) < E(IOO) . (28) 

Gluon loops corresponding to figure 1(b) yield a smaller contribution to the vacuum 
energy. To compute this energy we must use the gluon propagator suitably modified by the 
Meissner effect, which as we described above leads to the vacuum polarization P{kQ, k). We 
obtain 



Eg = I dkndk k^ In 



\ ^ Piko, k) 



kl + k^ 



(29) 



To estimate the result of this integral it is necessary to use the properties of P{kQ, k). Recall 
that P{ko, k) falls off like m|)Ao/fc at k » Aq, and similarly at large /cq [0- The dominant 
region of integration is therefore k^ ~ Aq and k^, << k << jj,, where A;* = m^ Aq is the 
momentum scale familiar from Landau damping. From this region of integration we obtain 

Eg ~ mlAlln{fi/k.) ~ gfi^Al , (30) 

since InA;,,, ~ 1/g. The result is parametrically smaller than the quark contribution. Note 
that this contribution to the energy is positive and hence prefers the least possible breaking 
of the color gauge symmetry. If this term were the dominant one then it would disfavor 
the formation of a condensate and the CFL vacuum would be the highest energy state! At 
asymptotic densities it is not the dominant term and the analysis from the quark loops 
stands. At lower densities, where the coupling is large, these contributions to the energy 
become more important but we lose control of the calculation. Their effect is to lower the 
energy gap between the CFL vacuum and disoriented states with different color and flavor 
symmetry breaking patterns. 
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The contribution we just computed to Eg still cannot differentiate between relative color 
rotations between the LL and RR diquark condensates. This is because the energy Eg 
depends on the sums of the squares of the gauge boson masses induced. In order to be 
sensitive to LR coupling effects, it is necessary that both LL and RR contributions to P{kQ, k) 
appear simultaneously in the vacuum energy contribution. The first such graph is that of 
figure 1(d), and it is of the form 



E^^ ~ / dkndk fc2 In 



^<m^^ 



(31) 



This integral is dominated by the region /cq ~ Aq, /;; ~ A;*, leading to the result 

This effect is further suppressed by a power of the coupling constant. We see that in the 
weak coupling limit the vacuum energy required to disorient the LL and RR condensates 
in color space is rather small. This suggests that even at asymptotic densities in the two 
flavor case it might be possible to disorient the diquark condensates from their lowest energy 
configuration. We have yet to determine what this lowest energy configuration is, and hence 
whether parity is violated in the two flavor case. In principle, one should minimize Eg^ as a 
function of the relative LL and RR color orientation. Instead, we will give a simple argument 
that the condensates prefer to ahgn. We noted in the last section that including the Meissner 
screening in the gluon propagator leads to a decrease in the gap size. This is a small effect 
at weak coupling, and was negligible compared to the color-flavor structure of the quark 
propagator. However, in determining LL-RR alignment it is the main effect. In the two flavor 
case none of the gluons responsible for the attractive interaction are Meissner screened, as 
long as the LL and RR condensates align. That is the quarks which condense are those that 
transform under the unbroken SU(2) subgroup of SU{3)c. However, any misalignment leads 
to the LL condensate screening the RR channel and vice versa, decreasing the condensates 
and thereby increasing the energy. Hence in the two flavor case the condensates prefer to 
align and parity is preserved. In the three flavor case CFL gives all of the gauge bosons a 
mass and this effect is absent. 

In both the two and three flavor cases there remains the possibility of parity violation 
through a phase associated with the U{1)a symmetry [Q, |^, |^. Only instanton effects (highly 
suppressed at asymptotic densities) can distinguish these vacua. At lower densities instanton 
effects are expected to strongly break the U{1)a symmetry, since the rj' mass is dominated 
by these effects at zero density. 
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4 Conclusions and Phenomenology 

In this paper we analyzed the possible ground states of QCD at asymptotic densities. We 
verified that in the two flavor case, the symmetry breaking pattern is SU(3)c — > SU(2)c, 
while in the three flavor case, color flavor locking has the lowest vacuum energy. In neither 
case is parity spontaneously violated until the density is strictly inflnite ||^, |^. 

Our analysis of the energy surface governing color superconductivity suggests possible 
experimental signatures in heavy ion collisions. In particular the existence of relatively 
flat directions along which color and flavor symmetry breaking patterns change raises the 
possibility of domains of disoriented condensates, each with distinct hadronization properties. 
In the two flavor case the LL and RR condensates, each of which break SU{3)c -^ SU{2)c, are 
only aligned by a subleading term in the vacuum energy calculation. In a heavy ion collision 
we might expect the condensates to be misaligned by an arbitrary SU{3)c transformation, 
leading to violation of parity and complete breaking of the color group. In the three flavor 
case we might expect much the same. Here the gauge loop contributions to the vacuum 
energy from gluon loops will tend to reduce the energy difference between the CFL and, 
for example, the (1,0,0) condensates as discussed above. The strange quark mass also tends 



to reduce the energy gap between these two condensates as discussed in |]T3[. For some 
(uncalculable) value of rris of order Aqcd we expect a phase transition between these two 
condensates as the number of light flavors changes from 3 to 2. Thus for realistic values 
of rris, and densities and temperatures achievable in heavy ion collisions, we might expect 
disoriented condensates to form with a range of possible color symmetry breaking patterns 
appearing on a collision by collision basis. 

To see how such variation in color symmetry breaking might be seen in an experiment 
we consider the extreme case where the SU{2)c subgroup is left unbroken (this is the true 
vacuum of the two flavor theory) . Consider a region which in the wake of a heavy ion collision 
volume is sufficiently cool and dense to allow the formation of a diquark condensate, with 
gauge symmetry broken to SU{2)c. The region presumably expands and cool in the usual 
fashion. However, one color of quark (e.g. red) does not participate in the condensation 
and its propagator is unaffected by A. It is also more weakly interacting since its color 
corresponds to precisely the broken part of the gauge group (gluons which couple to red 
quarks are screened by the Meissner effect). The remaining two colors of quarks participate 
in Cooper pairing and interact strongly with the plasma, so they do not disperse as quickly. 
The red quarks will therefore tend to flow to the surface of the flreball, providing a mechanism 
for macroscopic transport of color charge. Note that the condensate is stable under this 
charge separation since it is the condensate favored by an SU{2)c theory with two flavors. 
Furthermore, with SU{3)c broken to SU{2)c there is no restoring force which prevents this 
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charge polarization. On leaving the superconducting volume red quarks will suddenly be 
required to hadronize because their color charge can now support long range fields and they 
become aware of the large value of the other two color charges in the center. We expect that 
this color polarized fireball will hadronize very differently than a quark gluon plasma which 
is locally color neutral. Naively one expects quark anti-quark production on the boundary of 
the color charge separation in order to enforce charge neutrality. The separated red charge 
would then emerge as energetic hadrons, leaving a cooler central region behind. 

The scenario described above is the extreme case of a fully unbroken SU{2)c subgroup. On 
the other hand a CFL state treats all colors equally and there will be no charge polarization. 
On an event by event basis we expect variation between these two extremes. The most likely 
signal of such events is a departure from the standard thermal distribution so far observed 
in heavy ion collisions |2^, both on an event by event basis and averaged over many events. 



Acknowledgements: SH would like to thank D.K. Hong, V. Miransky, M. Nowak and 
M. Rho for useful discussions. Part of this work was performed at the Korean Institute 
for Advanced Study (KIAS). SH is supported under DOE contract DE-FG06-85ER40224. 
NJE acknowledges the support of a PPARC Advanced Fellowship. MS is supported under 
DOE contract DE-FG02-91ER40676. JH is supported in part by the Natural Sciences and 
Engineering Research Council of Canada and the Fonds pour la Formation de Chercheurs et 
I'Aide a la Recherche of Quebec. 



References 

[1] D. Bailin and A. Love, Nucl. Phys. B190 (1981) 175; Nucl. Phys. B190 (1981) 751; 
Nucl. Phys. B205 (1982) 119; Phys. Rep. 107 (1984) 325). 

[2] M. Alford, K. Rajagopal and F. Wilczek, Phys.Lett. B422, 247 (1998). 

[3] R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, Phys.Rev.Lett.81, 53 (1998). 

[4] R. Rapp, T. Schafer, E. V. Shuryak and M. Velkovsky, hep-ph/[A9904353. 

[5] N. Evans, S. D. H. Hsu and M. Schwetz, Nucl.Phys.B551, 275 (1999); Phys.Lett. B449, 
281 (1999). 



[6] T. Schafer and F. Wilczek, Phys.Lett.B450:325-331,1999; |hep-ph/9810509 
[7] D. T. Son, Phys.Rev.D54 094019 (1999). 



15 



[9] 
[10] 

[11] 
[12] 

[13] 

[14] 
[15] 

[16] 
[17] 
[18] 
[19] 
[20] 
[21] 



[22] 
[23] 
[24] 



R. D. Pisarski and D. H. Rischke, |nucl-th/981lT04 . 
R. D. Pisarski and D. H. Rischke, |nucl-th/ 9903023 . 
D. K. Hong, hep-pli/98125T0 



D. K. Hong, hep-ph/ 9905523 



D. K. Hong, M. Rho and I. Zahed, |hep-ph/990655T| . 

M. Alford, K. Rajagopal and F Wilczek, Nucl.Phys.B537, 443 (1999); M. Alford, J. 
Berges and K. Rajagopal, [hep-ph/9903502 . 



J. Berges and K. Rajagopal, Nucl.Phys.B538, 215 (1999). 



S. Hands, J.B. Ko gut, M.-P. Lomba rdo and S.E. Morrison, |hep-lat/ 990203^ ; S. Hands 
and S.E. Mo rrison |hcp-iat/99020l^ , hep-lat/990521. J.B. Kogut, M.A. Stephanov and 
D. Toublan, |hep-ph/9906346| ! 



T. Schafer and F. Wilczek, PRL 82 (1999) 3956; |hep-ph/99035nB| . 

D.K. Hong, V.A. Miransky, LA. Shovkovy and L.C.R. Wijewardhana, |hep-ph/9906478 

T. Schafer and F. Wilczek, |hep-ph/9906"5T2t 



R. D. Pisarski and D. H. Rischke, Qucl-th/9907041 



S.D.H. Hsu and M. Schwetz, |hep-ph/99083T0| . 



G. Benfatto and G. Gallavotti, J. Stat. Phys. 59, 541 (1990); Phys. Rev. C42 (1990) 
9967; R. Shankar, Physica A177, 530 (1991); Rev. Mod Phys. 66, 129 (1993); J. Polchin- 
ski, in Proceedings of the 1992 TASI, eds. J. Harvey and J. Polchinski (World Scientific, 
Singapore 1993). 

M. Le Bellac, Thermal Field Theory, Cambridge University Press, Cambridge, 1996. 

J.M. Cormwall, R. Jackiw and E. Tomboulis, PRDIO, 2428 (1974). 

F. Becattini, M. Gazdzicki, J. Sollfrank, Nucl. Phys. Proc. Suppl. 71 324 (1999); F. 
Becattini, M. Gazdzicki, J. Sollfrank, Nucl. Phys. A638 403 (1998). 



16 



